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Abstract. We report on our recent calculation of the inclusive direct photo- and hadroproduction of the J /iff meson at next- 
to-leading order within the factorization formalism of nonrelativistic QCD. We fit the color-octet (CO) long-distance matrix 

elements (ff J IV( x sf )), (ff J r /V( 3 sfl)) and (0 J /V( 3 P^ ] )} to the transverse momentum (p T ) distributions measured by CDF 
at Fermilab Tevatron and by HI at DESY HERA and show that they also successfully describe the pj distributions from 
PHENIX at BNL RHIC and CMS at the CERN LHC as well as the photon-proton cm. energy and (with worse agreement) 
the inelasticity distributions from HI. In all experiments, the CO processes are shown to be indispensable. 
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INTRODUCTION AND OVERVIEW 

The factorization formalism of nonrelativistic QCD 
(NRQCD) [1] provides a rigorous theoretical frame- 
work for the description of heavy-quarkonium produc- 
tion and decay. This implies a separation of process- 
dependent short-distance coefficients, to be calculated 
perturbatively as expansions in the strong-coupling con- 
stant 0£ s , from supposedly universal long-distance ma- 
trix elements (LDMEs), to be extracted from experiment. 
The relative importance of the latter can be estimated by 
means of velocity scaling rules; i.e., the LDMEs are pre- 
dicted to scale with a definite power of the heavy-quark 
(Q) velocity v in the limit v -C 1. In this way, the the- 
oretical predictions are organized as double expansions 
in 0Cj and v. A crucial feature of this formalism is that 
it takes into account the complete structure of the QQ 
Fock space, which is spanned by the states n — 2S + 1 ]jp 
with definite spin S, orbital angular momentum L, total 
angular momentum J, and color multiplicity a = 1 , 8. In 
particular, this formalism predicts the existence of color- 
octet (CO) processes in nature. This means that QQ pairs 
are produced at short distances in CO states and subse- 
quently evolve into physical, color-singlet (CS) quarko- 
nia by the nonperturbative emission of soft gluons. In the 
limit v — > 0, the traditional CS model (CSM) is recov- 
ered in the case of S- wave quarkonia. In the case of J/y/ 

production, the CSM prediction is based just on the 3 S^ 
CS state, while the leading relativistic corrections, of rel- 
ative order G(y A ), are built up by the l sf\ 3 sf\ and 3 pf ] 
(/ = 0,1,2) CO states. 
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The greatest success of NRQCD was that it was able 
to explain the J/y/ hadroproduction yield at the Fermi- 
lab Tevatron [2], while the CSM prediction lies orders of 
magnitudes below the data, even if the latter is evaluated 
at NLO [3, 4]. The situation is similar for the transverse 
momentum (pr) distribution at BNL RHIC [5]. Also in 
the case of J/y/ photoproduction at DESY HERA, the 
CSM cross section at NLO significantly falls short of 
the data [6, 7]. Complete NLO calculations for the CO 
contributions were performed for inclusive J/y/ produc- 
tion in two-photon collisions [8], e + e~ annihilation [9], 
and direct photoproduction [7]. As for hadroproduction 
at NLO, before this talk was given, the CO contributions 
due to intermediate and 3 sf^ states [4] were calcu- 
lated as well as the complete NLO corrections to %j pro- 
duction, including both 3 pj^ and 3 s[ 8 ' contributions [10]. 

In order to convincingly establish the CO mechanism 
and the LDME universality, it had been an urgent task 
to complete the NLO J/y/ hadroproduction calculation 
by including the full CO contributions. This was actually 
achieved in the work [11] presented at this conference. 

Our strategy for testing NRQCD factorization in J/y/ 
production at NLO is as follows. We first perform a com- 
mon fit of the CO LDMEs to the pj distributions mea- 
sured by CDF in hadroproduction at Tevatron Run II 
[12] and by HI in photoproduction at HERA1 [13] and 
HERA2 [14] (see Table 1 and Fig. 1). We then compare 
the p T distributions measured by PHENIX at RHIC [15] 
and CMS at the LHC [16] as well as the W and z dis- 
tributions measured by HI at HERA1 [13] and HERA2 
[14] with our respective NLO predictions based on these 
CO LDMEs (see Fig. 2). For details on the calculation 
and the input parameters used, we refer the reader to 
Ref. [11]. 



TABLE 1. NLO fit results for the J/ y CO LDMEs. 



4 ] )) 



(4.76 ±0.71) x 1(T 2 GeV 3 
(2.65 ±0.91) x 10~ 3 GeV 3 
(-1.32±0.35) x 10~ 2 GeV 5 



FIT TO HERA AND TEVATRON DATA 

The pr distribution of J/y hadroproduction measured 
experimentally flattens at pj < 3 GeV due to nonper- 
turbative effects, a feature that cannot be faithfully de- 
scribed by fixed-order perturbation theory. We, therefore, 
exclude the CDF data points with pj < 3 GeV from 
our fit. We have checked that our fit results depend only 
feebly on the precise location of this cutoff. The fit re- 
sults for the CO LDMEs corresponding to our default 
NLO NRQCD predictions are collected in Table 1. In 
Figs. 1(a) and (b), the latter (solid lines) are compared 
with the CDF [12] and HI [13, 14] data, respectively. 

For comparison, also the default predictions at LO 
(dashed lines) as well as those of the CSM at NLO 
(dot-dashed lines) and LO (dotted lines) are shown. In 
order to visualize the size of the NLO corrections to 
the hard-scattering cross sections, the LO predictions are 
evaluated with the same LDMEs. 

We observe from Fig. 1(c) that the 3 fj 8 ' short-distance 
cross section of hadroproduction receives sizable NLO 
corrections that even turn it negative at pj 7 GeV. 
This is, however, not problematic because a particular 
CO contribution represents an unphysical quantity de- 
pending on the NRQCD scale /x A and the choices of the 
renormalization scheme and is entitled to become nega- 
tive as long as the full cross section remains positive. 

In contrast to the situation at LO, the line shapes of 
the ^Sq and 3 pj 8 ' contributions significantly differ at 
NLO. Therefore we can now, in our combined HERA- 
Tevatron fit, independently determine [ff J l^ '(}$)) and 
(0 J /V( 3 pf ] )). Notice mat (tfJ/vfpW)) CO mes out neg- 
ative, which is not problematic for the same reasons as 
explained above for the short distance cross sections. 



PREDICTIONS FOR FURTHER DATA 
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FIGURE 1. NLO NRQCD predictions of J/y hadro- and 
photoproduction resulting from the fit compared to the CDF 
[12] and HI [13, 14] input data. 
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FIGURE 2. NLO NRQCD predictions of J/y hadro- and 
photoproduction resulting from the fit compared to RHIC [15], 
CMS [16], and HI [13, 14] data not included in the fit. 



We observe from Fig. 2 that our NLO NRQCD predic- 
tions nicely describe the pr distributions from PHENIX 
[15] (a) and CMS [16] (b) as well as the W distributions 
from HI [13, 14] (c), with most of the data points falling 
inside the yellow (light shaded) error band. The NLO 
NRQCD prediction of the z distribution (d) agrees with 
the HI data in the intermediate z range, but its slope ap- 
pears to be somewhat too steep at first sight. However, 
the contribution due to resolved photoproduction, which 



is not yet included here, is expected to fill the gap in the 
low-z range, precisely where the resolved contribution 
is peaked. Near the high-z endpoint region, the NRQCD 
expansion is understood to break down, and the NRQCD 
series could be resummed via the introduction of univer- 
sal shape functions [17], possibly in the context of soft 
collinear effective theory [18]. 



CONCLUSIONS 

We have performed a complete NLO calculation of the 
inclusive J/y direct photo- and hadroproduction cross 
sections including the full CS and CO contributions. 

We conclude from Figs. 1(a) and (b) and 2(a)-(d) 
that all experimental data sets considered here signifi- 
cantly overshoot the NLO CSM predictions, by many 
experimental standard deviations. Specifically, the ex- 
cess amounts to 1-2 orders of magnitude in the case of 
hadroproduction and typically a factor of 3 in the case 
of photoproduction. On the other hand, these data nicely 
agree with the NLO NRQCD predictions, apart from 
well-understood deviations in the case of the z distribu- 
tion of photoproduction. 

In our work for the first time a multi-process fit of the 
CO LDMEs is performed, which come out to be consis- 
tent with NRQCD scaling rules. Our work therefore con- 
stitutes the most rigorous evidence for the existence of 
CO processes in nature and the LDME universality since 
the introduction of the NRQCD factorization formalism 
15 years ago [1]. 

We should remark that our theoretical predictions refer 
to direct J/y production, while the CDF and CMS data 
include all prompt events and the HI and PHENIX data 
even non-prompt ones. However, the resulting error turns 
out to be small against our theoretical uncertainties and 
has no effect on our conclusions. 

We would like to thank Geoffrey Bodwin and Gustav 
Kramer for useful discussions and Barbara Jacak and Ce- 
sar Luiz da Silva for help with comparing our predictions 
to the PHENIX data [15]. 

COMPARISON WITH REE [19] 

As a last point, we compare our results with those ob- 
tained by the authors of Ref. [19], another full NLO 
NRQCD calculation of J / \j/ hadroproduction, which ap- 
peared after this talk was held. Adopting their inputs, we 

find agreement with their results for the 3 S^\ l sf \ 3 sf\ 

and 3 pJ 8 ' contributions. However, their fitting philosophy 
greatly differs from ours. Specifically, they only fit to 
Tevatron data with p T > 7 GeV, but account for prompt 
production, while we jointly fit to Tevatron data with 
Pt > 3 GeV and HERA data neglecting feed-down con- 
tributions. Detailed investigation reveals that the feed- 
down correction and the shift in the lower pj cut on 
the Tevatron data only moderately affect our joint fit. 
However, excluding the HERA data altogether renders 
the fit greatly underdetermined. Faced by this, Ma et al. 
[19] perform a constrained fit to just two linear combi- 
nations M and Mi of (ff J ^{ l sf)), (tflvQsf)), and 
(0 J /V( 3 P®)) defined in such a way that the 3 pf ] con- 



tribution is effectively accounted for by the l S ' and 3 S\ ' 
ones multiplied by Mo and Mi, respectively. Determining 
the very combinations Mo and M\ from a three-parameter 
fit just to the Tevatron data with pj>l GeV, we find that 
M\ has attached to it an error of almost 100%. This is be- 
cause Mo and M\ do not precisely correspond to eigen- 
vectors of the three-dimensional covariance matrix and 
the linear combination of (@ j lv( l sf )), (ff J ^{ 3 sf ] )), 

and (& J / V ( 3 P^)) corresponding to the third eigenvector 
carries a sizable error feeding into M\ . 
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